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A THEORETICAL/EXPERIMENTAL PROGRAM TO DEVELOP 
ACTIVE OPTICAL POLLUTION SENSORS 


By 

Frank S. Mills 1 and Roger N. Blais 2 
I. Introduction 

The intent of this research project was to develop and apply 
Light Detection and Ranging (LIDAR) technology to the assessment 
of air quality, and to evaluate its usefulness by actual field 
tests. Necessary hardware, to be described below (Section 3Aii) 
was successfully constructed and operated in the field. Measure- 
ments of necessary physical parameters, such as SO 2 absorption 
coefficients were successfully completed, and theoretical pre- 
dictions of differential absorption performance were reported. 
Plume modeling improvements were proposed (cf . appendix) . A 
full scale field test of equipment, data analysis and auxiliary 

\ 

data support was conducted in Maryland during September 1976. 
Thus, s i gnif icant_s trjLdes were made. in_all four areas (system 
development and demonstration; theoretical; field measurements; 
and modeling) described in the work statement of the original 
proposal. 

The following report will briefly summarize work previously 
reported, and then describe in detail the recent development in 
the work of Old Dominion University personnel both at Langley 
Research Center and on the Norfolk campus. 


1 Research Associate, Old Dominion University Research Foundation, 
Norfolk, Virginia 23508. 

2 Assistant Professor, Department of Physics and Geophysical 
Sciences, Old Dominion University, Norfolk, Virginia 23508. 



II . Work Previously Reported 


Material included in progress reports will be outlined 

here, but not reproduced. Three reports dated respectively' 

June 6, 1975; October, 1975; and December, 1975 were submitted. 

The June 6, 1975 report covered four basic topics 

reported in detailed appendices. First, a report on 

"Measurements of S0 2 Absorption Coefficients Using a Tunable 

Dye Laser" by R. T. Thompson, Jr., J. M‘. Hoell, Jr., and 

o. 

W. R.. Wade. Absorption coefficients for the 3001.8 A, 
o o 

2981.0 A and 29 62. A centered electronic-vibrational transitions 

o 

were reported with a. wavelength uncertainty of +0.1 A. Next, 
a report on- "Remote Sensing of Atmospheric SO 2 Using the 
Differential Absorption Lidar Technique" by Hoe-11, Wade and 
Thompson was included. Results from a computer simulation 
of a DIAL system indicated that commercially available 
technology could achieve measurement sensitivities less than 
2 ppb with spatial resolution of 500 m over ranges of less than 
2 km. Third, Thompson reported "Sensitivity Predictions for 
Differential Absorption and Scattering Lidar." General 
statistical error analysis equations were developed for 
analyzing the sensitivity of a DAS system, and they were 
evaluated for three experimental situations : ground level 

nitrogen dioxide (NO 2 ) measurements over a horizontal path; 
measurements of ozone (0^) depletion in jet engine wakes 
at 20 km altitude; and orbiting platform, nadir viewing 
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atmospheric ozone distribution measurements. Finally, 

Thompson and F. Allario reported on "Optical Properties of 
Tunable Diode Laser Radiation," in which divergence, 
polarization and pattern anomalies in the far field pattern 
of tunable dye lasers were investigated. 

The October, 1975 report (Technical Report PGSTR-PH75-12) 
by S. K. Poultney, M. L. Brumfield and J. S. Siviter was 
subtitled Quantitative Remote Raman Lidar Measurements of 
Pollutants from Stationary Sources . A detailed study of 
using Raman lidar was conducted using a calibration tank at 
LaRC. It was shown that typical stack exit- concentrations 
of 500 ppm S0 2 could be measured to an accuracy of 10 percent 
at a distance of 300 m, with integration times of 30 minutes. 

The December, 1975 report by R.- T. Thompson, Jr., 
"Differential Absorption and Scattering Sensitivity Predictions" 
was a much fuller account of the studies summarized in the 


June 6 report. 



Ill . Work Previously Unreported 

Included below are results previously unreported. 

Section A consists of work done by ODU personnel at LaRC, 
and it is divided into two subsections, the first theoretical 
and the second experimental. It is the work of Dr. Mills 
and. his group. Section B contains a brief description 
of various activities of ODU personnel working on campus. 

A fuller presentation of their results is found in the 
appendix. 

A. Work Done At Langley Research Center 

i) Theoretical Support for the. NASA Langley Water Vapor 
DIAL Experiment. 

The purpose of the experiment was to demonstrate the 
feasibility of the differential absorption lidar technique 
for measuring water vapor profiles in the troposphere up to 
3 km. For optimum operation the technique should use water 
vapor lines which are relatively isolated in frequency, have 
an absorption cross-section which is relatively insensitive 
to changes in temperature, and have an absorption cross- 
section which is large enough to produce measurable 
absorption, but not large enough to completely absorb the 
laser radiation. The. work described here involves determining 
what water vapor lines would be least sensitive to temperature 
changes . 


REPRODUCIBILITY OF THE 
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Assuming an isolated Lorentz shaped absorption line, the 
following expression for the absorption coefficient, k- 1 , may 
be written. 



S Y 

ir(A\>2 + 


( 1 ) 


where S is the line strength, y is the Lorentz half-width, 
and Av is the difference between the frequency of interest 
and the frequency of the absorption line. Both -the line 
strength S and Lorentz half-width y are temperature dependent. 
The half-width temperature dependence can be expressed by 


the 


following expression. 



( 2 ) 


where y q is the defined as the half-width at temperature T q 
and T is the temperature. The temperature dependence of the 
line strength can be expressed as follows': ' • 

S = Kq T-3/2 e -hcE/kT 


where Kq is a constant, T is the temperature, h is Planck's 
constant, c is the speed of light, k is Boltzmann’s constant, 
and E is the energy of the lower energy level of the transition 
expressed as wave number. 

The dye laser used in the DIAL experiment has a line width 
which is comparable to the absorption line width. Therefore, 
the absorption which is measured by the DIAL is the absorption 
integrated over the laser line width. The approach used 
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here will be to calculate the integrated absorption assuming 
some laser line shape, and then determine for a given 
temperature the lower state energy for which the integrated 
absorption is least sensitive to temperature variations. 

Calculations were made for two different laser line shapes, 
square and triangular. The actual laser line shape is 
somewhere in between. Let F(v') be a function describing 
the laser line shape where 

co 

^ . / F(v’) dv ' = 1 (4) 

' —00 

Figure 1 shows the two laser line shapes graphically. 

For the square line shape F (v * ) may be written 


F (v 1 ) ’ 

= l/2w 

-w < v ’ < w 

(5) 

F(v ') 

= 0 

1 V 1 I > w 



For the triangular line shape F(v') may be written 

F(v') = l/2w - v*/4w^ 0 £ v ' < 2w (6) 

F(v’) = l/2w + v'/4w2 -2w < v’ <_ 0 

F ( v ' ) = 0 jv ' J > 2w 

For the square laser line shape, the integrated absorption 
k' is 

=-2/ S y (l/2w ) dv ' ( 7 ) 

0 it (v 1 ^ + y^) 

or 


k' = (S/wir) tan ^(w/y) 


( 8 ) 



For the triangular laser line shape , the integrated 


absorption k 1 is 


k , „ 2/ 2W S y (l/2w ~ v * / 4w 2 ) dv 
0 x (v '. 2 + y 2 ) 


( 9 ) 


or 


k' - (Sy/irw) jcj./Y)tan“l (2w/y) -Cl/4w)ln (1 + 4w 2 /y 2 7J (10) 


Now , for each laser line shape the temperature dependence 
of S and y from equations 2 and 3 can be substituted into 
the expressions for the integrated absorption, and the lower 
state energy for temperature insensitive lines can be 
calculated by taking the derivative of the integrated 
absorption with respect to temperature, setting the 
derivative equal to zero and solving for E. 

Following the procedure described above for a square 
laser line, the integrated absorption- k ' in equation 8 can 
be written 


k , . M 3/2 e - hcE/kT tan - 1 ( . 

wit Y 0 'V 

or 

k' = K , T _3 / 2 e" hcE / kT tan" 1 {wTVYo T o J5 5 

where 


( 11 ) 


( 12 ) 


K' = K 0 /wit . 


(13) 
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Taking the derivative of equation (12) with respect to .T 


*1 „ K . T -2 e -hoEAI *121. T'^tan 1 ( WT , ) 
dT (_kT y QT h 


3 T -l5 tan -l ( _»?^) 

W 


2y 0 T 0 = ( 1 + Jd|_) 

V T o 


For temperature insensitive lines = 0. 


hS* T^tan-l ( wT 2 ' ) _3 T -^ tan ~l ( wT 3 } + w . » 

v T ^ 2 v m ^2 2vT^i 2 

Y 0 T 0 y 0 T 0 2y 0 T 0 .1+ w^T . 

( 2„ J 

Y 0 T 0 


tan-1 (-SSI.) = 3 tah -l ( _w** , _ 


WT^- 


Y 0 T 0 


V 2 


H 2y 0 T 0 % i l+ w 2 T j 

v 0 2t o 


Solving for E 


« _ kT 3 _ 
E " h? 2 


2 Y .T% tan" 1 (_wf!_) 

V 0 2t 0 YoTo* 5 


Following the. same procedures for the triangular laser 
line, the integrated absorption k* in equation (10) can be. 


k , = KqToTq^T- 2 e -hcE/kT 1T^ 

TTW Y T " 

o o 


tan (2 wT^_ ) 

■*o T o h 


- JL In (1 + T ) 

4w v 2m 

Y o r o 
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or 


k' = K'T 


, q ,-3/2 e -h=E/kT tan -l ( WL ) 




0 X 0 


- *L Y0 T 0 ^2 e ^ hcE A%(i + 


■y 2 T 
0 0 


) 


.where 


K 


K' = 


_0 

TTW 


(19) 

(20) 


Taking the derivative of equation 19 with respect to T 


||L . K ,^-hcEA4-5/2 


% 

hcE . _ i / 2wT 2 > 3 . -ic 2wT \ 

Et" tan (— TV - 2 tan 1 

t 0 V 


- *0*0 


wT' 


Y 0 T 0 !g (X+ '*”^ . . ) 

V T 0 


' - K^ v T ^ e -hoE/kT T -3 


r 




\ _ 


hcS ln(l + Jg 2 T) - 2 In (1 + 4.y! 2 ? .) 

2m 

Y o T o 


Y 2 T 
T 0 0 


4w 2 T 


Y 2t„ (1 + 4 ^) 


'0 0 


dk 1 


y 0 2t 0 


For temperature insensitive lines = 0 


kT 


tan (- 


Y 0 T 0 


h 

^ “” tan ( 


hcE J.-1 / 2wT^v 3 , ~1 , 2WT 35 


hcE ^0 H , ,, , 4w 2 T v Y 0 T 0 , 4w 2 T y 

vt iT" ln(l+ — 5 — ) T ln(l+ 2m } 

kT 4wir^ Y 0 2 T 0 2wT"* Y q T 


Y 0 T 0 


( 21 ) 


( 22 ) 


or 
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hcE 

kT 


tan^^L)- I0^ ln(1+ 4A } 


Y o T o' 


4wT i 2 


Y 0 - T 0 


= itan-l (2 "tL) - Y ° T ° ln d+i^ 


vo T o 


2wT^ 


Y 

Y 0 *0 


) 


( 23 ) 


or 


E = 


kT 

fic 


3 

2 


« in a + -ia!s_ 

8wT Y 0 T 0 


tan-1 (2wTl) - 2^1 In (It 42Z?*-) 


y 0 T 0 


4wT : 


2 m 
Y 0 T 0 


(24) 


J 


A sample calculation will be made for each laser line shape 
assuming that w, the laser line half-width and y^, the absorption 
line half -width, are equal, and that T = Tq = 296K. For the 
square line shape using equation 17 , the result is 


E = 1.18 k'E = 243 cm" 1 
he 


(25) 


For the triangular line shape using equation 24, the result is 


E = 1.21 S a 249 cm" 1 
he 


(26) 


The results’ of the sample calculations indicate that if 
the laser line width is comparable to the absorption line width 
the lower state energy for temperature insensitive lines is 
relatively independent of the exact laser line shape. 

Using a procedure similar to that outlined above, it can 
be shown that for a monochromatic laser line , the lower state 
energy for temperature insensitive absorption lines is just 
kT/hc or for T = 39 6K, 206 cm" 1 . Thus, the fact that the 
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laser line has a finite width has more effect on determining 
the lower state energy for temperature insensitive lines than 
the exact shape of the laser line. 

ii) Analysis of the DIAL Technique for Remote Probing of 
SO 2 in the Atmosphere 

SO2 from fossil, fuel powerplants is one of the primary 

problems contributing to degradation of regional air quality 

in this country. It is, therefore, desirable to develop 

remote sensing techniques which can monitor SO2 emissions 

from stationary sources or ambient SO2 concentration in an 

1 

urban environment. Poultney, et al. , have reported the 
results of an experimental program to determine the usefulness 
of the Raman lidar as a technique for remotely monitoring 
SO2 emissions from fossil fuel powerplants. In this report, 
currently available S 0 2 spectra are used to predict the 
performance of the Differential Absorption Lidar (DIAL) 
technique operating near 4 ym and 300 nm for remote sensing 
of SC>2 in the atmosphere. These two wavelength regions are 
covered by lasers which NASA Langley Research Center currently 
has under development. 

The DIAL concept can be described qualitatively as 
follows. Pulsed laser radiation at two wavelengths is 
transmitted into the atmosphere. The two wavelengths are 
selected so that one, called the on wavelength, is absorbed 
by the gas of interest and the other, called the off wave- 
length, is not. The backscatter return signal as a function 
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of range at each wavelength is collected by an optical 
receiver. Now define a range cell as the distance from 
% to R2. The average concentration of the absorbing gas 
in the range cell may be determined from the on and off 
wavelength returns at R-j_ and ■ The situation is 
analagous to a dual -beam spectroscopy experiment where the 
off wavelength return corresponds to the reference beam, ' 
the on wavelength corresponds to the sample, beam, and the 
range cell corresponds to the sample cell. The transmittance 
of the absorbing gas in the range cell is then just the 
ratio of the on wavelength returns at R2 and R^ divided by 
the ratio of the offv wave length returns at R2 and R^ • 

The quantitative. expression for absorber concentration 
can be found using the lidar equation: 


P r (R,l,t) = k ^Lp 0 (X) 6 (R, X , t) 


x exp {- 2 /g |jx sc (r,X,t) + N A {r,t) a A (X,P,T) 

+ N int (r ' t)a int(X,P,T)^] dr}. (27) 

Here P r is the received power as a function of range R, 
wavelength X, and time t; K is a constant; L is the laser 
pulse length; A is the receiver area; P Q {X), the transmitted 
power, is a function of wavelength; g is the backs catter 
cross section as a function of range, wavelength, and time; 
a Sc is the extinction coefficient for scattering as a function 
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of range , wavelength and time; N A is- the concentration of 
the absorbing gas of interest as a function of range and 
time; a A is- the absorption cross section of the gas of 
interest as a function of wavelength, pressure, P, and 
temperature, T; is concentration of any interfering 

absorbing gas; and u^nt is the absorption cross section of 
the interfering gas. 

A number of assumptions are made to simplify the analysis. 
The scattering terms are assumed equal for the on wavelength 
and the off wavelength returns, that is 6 (R,Aoff,t off) = 
g(R,Xon,t on) and a sc (R,Xoff, t off) -= a gc (R,Aon,t on). 

Also, the concentrations of the absorbing and interfering 
gases are assumed to remain constant between the off 
wavelength and on wavelength returns, that is, N A (R,t off) = 
N A (R,t on) and N^ n -t{R,t off) = Ni nt {R,t on) . The on and off 
wavelength absorption cross section for the gas of interest, 
c A (Aoff ,P ,T) and cr A ( Xon,P,T) , are assumed to be known and 
constant between R^ and R2. If (Aoff ,P,T) is not equal 

to ( Aon,P ,T) , then both the concentration and absorption 

cross section of the interfering gas must be known or must 
be determined by a separate experiment. 

With the above assumptions the average concentration 
N A of the gas of interest in the range cell between R^ 
and R2 is 

In 


N A = 


2(R2~R]_) £o A (Aon) -c A (Xoff)~ 


Pr on( R l> x Pr off( R 2 \ 
Pr off (%) X Pr on {R 2 ) 


( 28 ) 
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The quantities in equation 2 8 which are measured 
experimentally. . are the on and off. wavelength returns at R^_ 
and R 2 . Equation 28 can be rewritten as follows. 

Pr on(R2)/Pr Qn (Rl) _ £-2 (R 9 -R 7 )N a |g a (Xon) -a a (loff )]> (29) 

The term on the left is just the transmittance, of the gas of 
interest in the range cell from to R 2 . The exponent is 
the optical depth of the absorbing gas in the range cell. 

So from the. error in . determining the transmittance , which 
is the experimentally measured quantity, the corresponding 
error in determining .the optical depth, can be determined, and, 
assuming that the absorption cross section is known, the 
error in determining the concentration. 

The procedure then will be to select appropriate on and 
off wavelengths using available SO 2 spectra in the. 4 jim and 
300 nm wavelength regions. The absorption cross sections for 
the on and off wavelength will be used to make a plot of 
error in optical depth or absorber concentration as a function 
of optical depth or absorber thickness for a given error in 
measuring the transmittance where the absorber thickness is 

2 (r 2 -r 1 )n a . 

The spectrum of the v-j_+V 3 combination band of S0 2 near 

2 3 4 

4 ym has been studied by a number of workers ' ' using 

i 4 

grating spectrometers with resolutions ranging from . 017 cm J - 

-t 2 5 

to .48 cm - - 1 - . Recently, Pine has studied the band 

using a cw difference-frequency spectrometer with a resolution 
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of .0003 cm Figure 2 shows an absorption coefficient 
spectrum, observed by Pine , under atmospheric conditions. 

Using figure 2, a line pair of 2498.5 cm“^ and 2500.8 cm“^ 
was chosen for analysis, with the on wavelength line having 
a frequency of 2498.5 cm - ^ and the off wavelength line having 
a -frequency of 2500.8 cirT'*'.. The on wavelength absorption 
cross section is .5 (m% atm) and the off wavelength 
absorption cross section is 0 within Pine's error of 
+ .02 (m% atm)“-*-. 

g 

The spectrum of SO 2 near 300 nm was measured by Thompson 
with a resolution of ..02- nm. Figure 3 shows an absorption 
coefficient spectrum observed by Thompson 0 . Using figure 3, 
a line pair of 299.4 nm and 300.0 nm was chosen. The on 
wavelength, 300.0 nm, has an absorption cross section of 
33.1 (atm-cm)"*^. The. off wavelength, 299.4 nm, has an 
absorption cross section of 7.28 (atm-cm) The differential 
absorption cross section Qr A (A on ) - ^A^off^J t ^ ien 
25.8 (atm-cm) 

Figure 4 shows error in optical depth or absorber 
concentration as a function of optical depth for errors in 
measuring the transmittance of 1, 2, 5, and 10 percent. Also 
shown on the abscissa are the absorber thicknesses corresponding 
to the 4 ym differential absorption, the 300 nm differential 
absorption, and the 299.4 nm off line absorption. 
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Note that for any given error in measuring the trans- 
mittance/ the minimum error in the concentration determined 
from that measurement occurs for that combination of 
concentration and range cell for which the optical depth 
is equal to 1. Also the minimum error in determined 
concentration is approximately three times the error in 
measuring the transmittance. 

Assuming a range cell of 100 meters and a 5% error in 
measuring the transmittance , at 4 pm the minimum detectable 
average concentration in the range cell is about 5% parts 
per million and apparently the maximum detectable concentration 
is approximately 300 parts per million. The upper limit is 
actually much higher since the on wavelength could be tuned 
to a point where the absorption coefficient and therefore 
the optical depth are lower. 

At 300 nm, again assuming a 100-meter range cell' and~5%" 
error in measuring transmittance, the minimum detectable 
concentration is about 110 parts per billion. The maximum 
detectable concentration is determined, by the point where 
the absorption of the off wavelength line is so great that the 
return from the' far side of the range cell cannot be measured 
accurately. For a 100-meter range cell and 5% error in 
measuring transmittance, the maximum detectable concentration 
is 20 ppm. Thus, for measurements of powerplant emissions, 
the DIAL instruments operating at 300 nm and 4 ym are 
complementary . 
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Assuming 5% error in measuring transmittance, a one 
kilometer- column content measurement of SC >2 concentration 
as low as 11 parts per billion could be made using the 
DIAL technique at 300 nm. 

It should be emphasized that this analysis has not 
considered what accuracy in measuring transmittance is 
achievable. For any specific system, this would have to 
be determined from a consideration of specific system 
parameters such as that made by Thompson 7 . 
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Experimental Support for NASA Langley Plume Dispersion Program 

In the early part of 1976 an agreement was. made by 
NASA Langley Research Center and the Maryland Power Plant 
Siting Program to perform a joint experiment to determine 
the utility of the lidar as a method fox characterizing 
plume rise and plume dispersion from power plants. In - 
support of that experiment, a mobile lidar system was built 
using NASA equipment by personnel from Old Dominion University, 
NASA, and Wy-le Laboratories. 

Operation of the system was checked during late July 
by observing the plume from .the VEPCO plant at Yorktown, • • 
Virginia. 

In this report, the lidar system will be described. 

A later -report will describe the design of the joint field 
experiment and the results of that experiment. 

The lidar system consists of a ruby laser, a telescope 
receiver, a detector package, and associated instrumentation. 
The laser, telescope, and detector package are mounted on a 
searchlight-type mount with tracking capability. The entire 
system, including instrumentation and the searchlight mount, 
is contained on a flat-bed trailer. 

The laser used in the system is a Holobeam 600 
Q-switched ruby laser with a beam divergence of 3% milli- 
radians. The beam divergence was reduced to approximately 
1 milliradian by using an up-collimating telescope at the 
output of the laser. The output energy of the laser can 
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vary from .75 to 2.0 joules with a pulse length of 30 
nanoseconds. Thus, the ultimate range resolution attainable 
with the system is 4.5 meters. 

The receiver is a 12-inch Cassegrain type telescope 
with a field of view of approximately 4 milliradians . The 
detector package has provision for mounting two photo- 
multiplier tubes for extended dynamic range . For the plume 
dispersion experiment, only one tube was used. An RCA 
7265 photomultiplier tube was selected for use in the system 
since it was sensitive to the 694 nanometer laser radiation, 
and it could be easily gated to prevent overload from the 
close-in return . 

The instrumentation used consists of a high' voltage 
power supply and gating circuit for the photomultiplier 
tube, pulse generators used for timing, and the data 
acquisition system. 

The data acquisition system is based on a Digital 
Equipment Corporation model PDPll/10 minicomputer. The 
lidar return signal from the photomultiplier is recorded 
by a Biomation 8100 transient digitizer which can record 
2000 8-bit words of data at sample rates up to 100 MH Z . 

The Biomation has an analog output so that the data can be 
displayed continuously on an oscilloscope and a digital 
output from which the computer accepts the recorded data. 
Other data accepted by the computer include laser energy. 
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laser shot counter reading, and elevation and azimuth angle 
from the searchlight mount. The data from each laser firing 
is immediately recorded on magnetic tape for later processing 
using either the PDP 11/10 or the main Langley Research 
Center computer facility. Also attached to the minicomputer' 
is a Ramtek Graphics display which can display information 
in 16 shades of gray on an ordinary black and white 
television monitor. 

The computer has been programmed to allow data 
processing to proceed simultaneously with data recording. - 
This permits preliminary data analysis (such as range 
correction), and display on. a nearby real-time basis. 

Three different types of data display are available. One 
is called an A-scope display and is simply an x-y display 
with no intensity modulation of the display. Another type 
of display, called a z-scope display, is useful for 
applications where the lidar is pointed vertically. This 
display plots intensity versus height on a vertical line 
using the 16-shade gray scale modulation. The third type 
of display is the RHI or range, height, intensity display 
which is used for displaying plume dispersion lidar returns. 
For this display, the horizontal axis corresponds to the 
horizontal distance from the lidar and the vertical axis 
corresponds to height above the lidar. Intensity is then 
plotted along a line which corresponds to the lidar 
elevation angle. 
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The lidar was tested during late July by observing the 
plume from the VEPCO plant at Yorktown. The plant is oil 
fired with no precipitators on the exhaust. The stack is 
550 feet high. The lidar was located approximately 5.5 km 
from the stack with a viewing window about 30° to either 
side . The plume- was easily observable at a range of 8 km 
from the lidar and at a distance of 4 km downwind from the 
plant. Because of the restricted viewing window, the 
maximum useful range of the lidar could not be determined, 
but it is certainly greater than 8 km. 

During the test, malfunctions in some of the computer 
peripherals occurred and shortcomings in the software were 
found. These problems were all corrected before the joint 
experiment in Maryland i 



B. Work Done by ODU Resident Faculty. 

ODU faculty- devoted themselves to three main support 
efforts. First, Dr. Kindle and Dr. Blais consulted with 
Langley personnel on feasibility of field test procedures. 
Second, information was assembled on local climatology and 
on stack parameters of local sources for support of the 
Yorktown plume field test. Third, Dr. Blais initiated 
development of an adaptation to the Gaussian plume model 
using a surface interaction parameter first suggested by 
Csanady*. Geometric characteristics capable of lidar 
detection were proposed. Work is continuing on an error 
and sensitivity analysis of the modified model to make it 
more suitable for field application. Dr. Blais’, report 
constitutes the appendix. 


*G. T. Csanady, Aust. J. Phys. 8_, 545-550 (1955) . 
G. T. Csanady, Aust. J. Phys. 10, 559-564 (1957). 
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Figure 1. Laser line shapes. 
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APPENDIX 


LIDAR DETECTABLE GEOMETRIC CHARACTERISTICS 
OF GAUSSIAN PLUME MODELS WITH VARIABLE 
SURFACE BOUNCE 

' - ’ . Roger- N. Blais 

Department of Physics and Geophysical Sciences 
Old Dominion University . 

Norfolk, Virginia 23508, U.S.A. 

Abstract - A Gaussian Plume Modal is proposed containing 
a dimensionless bounce parameter, 3/ the fraction of 
material dispersed to the earth's surface that is reflected 
into the plume. The 6 = 0 case yields the free space 
plume in the z > 0 domain, with perfect trapping of effluent 
on the z .= 0 plane. The 6 = 1 case reproduces the Pasquill- 
Gifford Model with all effluent reflected a t t he surface . 

The .model is formulated to allow estimation of 6 from lidar 
measurements of plumes over various textures of flat 
terrain (e.g. forest, grass) . Lidar detectable geometric 
characteristics of plumes are described as functions of a 
non-dimensional downwind distance parameter, a 2 /H, and of 
6 as it ranges from 0 to 1. These characteristics are the 
altitude- of maximum concentration, and the mean altitude of 
the concentration . The measured mean altitude * s dependence 
on the maximum - altitude of lidar sampling is discussed. The 
influence of 0 on the total downwind flux of effluent and on 
surface concentration is described. 



NOMENCLATURE 


•C(x,y,z) concentration of effluent, kg m' 


H _ effective emission height, the sum of stack 

height and plume rise, m 

M n (y) n tia moment of n integrated over Aq from n = 0 

to r => u, see Equation 3,. dimensionless 


Q 


Qi 

u 


x,y ,z 


S 

n ' 

^m 

% 

A(a,y,n) 

Ao(dfn) 

y 

a y (x) 

cj z (x) 


rate of effluent emission at source (0,0,H), 
kg s” 1 


rate of effluent emission by virtual 
at (0,0, -H) , kg s”-*- 


imace source 


mean- wind speed, 


m s 


-1 


downwind, crosswind and vertical coordinates 
respectively, m 


normalized vertical plume thickness, dimensionless 

the dimensionless surface bounce parameter such 
that 0 <• 3 <_ 1. * 

normalized altitude, dimensionless 

t) at which maximum A Q occurs, dimensionless 

mean (expectation) value of n averaged over Aq 
from n = 0 to n = u 

normalized concentration of effluent, see Equation 
4, dimensionless 

normalized concentration of effluent on y = 0 
plane, see Equation 5, dimensionless 

the maximum normalized altitude to which lidar 
measurements are made, dimensionless 

standard deviation of'C(x,y,z) relative to y, m 

standard deviation of C(x,y,z) relative r to z in 
free space, m 


$(ct,8 ) 


total normalized mass flux through a vertical plan 
normal to the wind, see Equation 16, dijnensionless 



LIST OF CAPTIONS 


Fig. 1. the altitude of maximum Ag, as a function of 

normalized downwind distance parameter, a. 

Fig. 2. A.ltitude of the true mean (j\ m ) , of a measured mean 
with sampling interval 0 _< ri <_ 3 (rj^) , and the . 
altitude of maximum normalized concentration (n rc ) . 
Fig. 3. The normalized flux ($(a)) or fraction of original 
effluent mass still airborne at a, and the 

00 

normalization factor (Mg (°°) ) . 
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INTRODUCTION • 

• The 'growing -use of the Gaussian Plume Model (Csanady, 

1973; Seinfeld/ 1975; Sutton, 1932; Turner, 1970) for the 
design- of stacks, for the preparation- of environmental 
impact statements, and for the development of multiple 
source models (TRW, 1969) makes the validity of the 
assumptions upon which the model rests a legitimate economic 
and public health concern. Chamberlain (1966) , Bessemoulin 
.(1974) , Heines (1974) , Ragland (1975) , and Seguin (1973) , 
among others , have examined the boundary conditions at the 
earth's surface! The influence of the surface on the plume 
is amenable to lid-ar (laser radar) study if a simple field 
model can be used to relate vertical concentration profiles 
to a parameter that describes the surface. Lidar is better 
able to establish plume geometric characteristics than to 
define absolute mass concentrations, due to the difficulties 
of the Mie scattering problem for polydisperse aerosols . 

Field workers are seldom able to measure all the necessary 
input variables, which include the complex index of refraction 
for each of the various constituents, and their particle size • 
and shape distributions (Diermendjian, 1969) Geometric' 

characteristics of a plume, however ,- are readily defined 

\ 

under simple assumptions of particle homogeneity and isotropy 
in the spatial domain. 
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In addition , lidar systems can time average data over 
periods consistent with the assumptions of the Gaussian 
Models-. In this regard, lidar is superior to instantaneous 
photogrammetric techniques (Blais, 1975). 

This paper considers lidar detectable plume geometric 
characteristics capable of illuminating surface interactions, 
by exploiting a simple generalization of the Pascuill-Gifford 
Model (Gifford, 1961; Hay, 1957; Pasquill, 1961). As in the 
standard model, the surface is required to be relatively 
flat, but its texture may vary widely in effective porosity. 

The generalization is based on a surface bounce parameter. 
Non-dimensional expressions for altitude of maximum concentration, 
for true and for measured mean altitude of concentration, and 
for downwind flux are presented as functions of a nOn- 
dimensional downwind distance parameter, and of the surface 
bounce parameter. The influence of 8 on surface level 
concentrations is discussed. 
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THE BASIC MODEL' 

The Pasquill-Gif ford Model is modified by including 
a dimensionless bounce parameter, f3 , which ranges between 
zero and one. It multiplies the surface reflectance term 
such that: 

C = ( Q ) exp (-. Y ) [exp — ) +.6 exp (--k±SL_)l {1) 

2irUayC? z 2Cy 2 2o z 2 2c z 2 • 

In the case 3=1, presumed to be typical of a smooth flat 
surface, all of the material dispersed downward to the 
earth's surface remains airborne because it is bounced or 
reflected back into the plume. This is the condition for 
the regular Pasquill-Gifford model. The opposite extreme 
case, 3=0, presumed to be. approximated by dense forest, 
implies total absorption of material at the earth’s surface. 
The plume then behaves as a free space plume in the domain 
z > 0. 

The S parameter arises from the artifice of using a 
virtual image source located .at .x = y = 0 , z = -H to account 
for reflectance at the earth's surface. Usually, the rate 
of image source emission, Qi, is assumed to be equal to 
the real source emission rate, Q. Because the surface 
concentration depends upon summing the effects of the real 

r 

and the image source, and because Qi > Q is not physically 
meaningful, the Qi = Q case represents the highest possible 
surface concentration at any given field point. This worst 
case is useful in predicting the environmental impact of 



proposed smokestacks. The true influence of the surface 
'■ on plume material is undoubtedly between the extremes 
represented by 3 = 1- and B = 0, where/ mathematically 8 is 
defined as: ■ • 

S - Q ± /Q, • (2) 

though physically it. is dependent on surface properties. 

The presumed illustrative examples of~ "smooth flat .land 
implying 3-1/ and of dense forest implying- 3 — -0, bear 
further comment. For a single value of 3 to characterize 
plume dispersal requires that the terrain J_a. flat, and 
uniform, and that the plume is not influenced by complicating 
factors like penetration of the inversion, or trapping by - 
the top of the mixing layer. Given such conditions , two 
surfaces may be compared, both of which are flat and smooth,- • 
but which differ in porosity. The first surface is exemplified 
by a smooth hard surface, like pavement, or~ closely mowed 
grass. The second surface is dense forest canopy. The forest 
has a stagnant air layer beneath the canopy, and a comparative! 
enlarged ratio of absorbing surface area to land surface area. 
Thus,, it is anticipated that forest will be a poorer reflector, 
and a better absorber of effluent than pavement, or' mowed 
grass. Consequently, one expects 8 forest ..< 8 pav amen t - 
Lidar studies can determine values for 3 over diverse , but 

t 

uncomplicated, surfaces like water, swamp or moorland. 



( 5 ) 


THE NORMALIZED MODEL 

It is useful to formulate the model non-dimensionally , 
so that- it is - independent of particular functional forms 
of Cy(x) or o z (x) . The following dimensionless quantities 
are defined. 



= z_ As 2 ttUCvCH 
H ' Q 


( 3 ) 


Note that a z generally increases with x (Turner, 1970) and 
H is constant in the region of interest far enough downwind 

o 

for surface bounce to be significant. Consequently, a is 
a measure of x, in that ct 0 as x -► 0 , and a -> « mono.tonically 
as x -*■<». The relationship is non-linear, however, and 
depends upon stability class. 

Note also that A is slightly different from the usual 
normalized concentration parameter, CU/Q, which is dimensionally 
a reciprocal area. As will be seen, A is more convenient 
in the present formulation, because it absorbs the a y 
dependence of the model. 

Substituting Equations 3 into Equation 1, and isolating . 

3 by the extraction of a common monomial factor yields 


A(a,y,n) = a 1 exp(-_zi - [ exp ( 2n ) + b] . (4) 

2a y 2 2a^ a 2 

Next, defining A 0 (oc,n) = A(a,0,n) yields » 

A 0 (a,n) = exp,.. [n-t-lf ) f exp (2n) + (5) 

2a 2 a 2 

Thus, Aq represents the non-dimensional concentration on the 
plume vertical plane of symmetry. There are three reasons 
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for dealing, with Aq rather than with ... First, a lidar 
system samples along an entire line of sight. • In the 
downwind region where the plume is sufficiently tenuous 
for multiple optical scattering events to be insignificant, 
but for surface bounce to be influential, one can find 
the centroid of the plume and plot vertical profiles 
through that plane only. Next, even if one wishes to 
examine the total optical brightness of the plume due to 
diffuse scattered light, an integration of total mass- 
load. A, along an infinite horizontal cross-wind line of 
sight (from y =' to y = +°°) would differ from Aq only 
by the multiplicative factor /2tT G y. Finally, some .regional 
multiple source models, like the AQDM (TRW, 1969-)-, use a 
Gaussian vertical dependence, consistent with Aq,. but a 
non-Gaussian horizontal dependence incompatible with A. 
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GEOMETRIC PROPERTIES OF A g 

To use the above model for estimating 0 , the lidar 
system in the field must acquire the- following data. 

Profiles of Ag versus n are determined at a given value of 
a. • To accomplish this, the lidar must aim in a horizontal 
cross-wind direction, and scan vertically between n = 0 
and a maximum normalized altitude r) - y. Acquired data 
must be stored and the lidar then changes azimuth to make 
a' 'second Ag versus n profile at a new a. After measuring 
profiles for several a values , the entire sequence is * 
repeated, until time averaged profiles are established for 
each a. Averaging times should be commensurate with the 
assumptions underlying the Gaussian Model (Slade, 1968) . . 

Two physical properties of the profiles are then established 
r| m , the altitude of maximum Ag? and rfy , the mean altitude 
or expectation value of ti over the- Ag(ct,n) distribution 
measured from n = 0tof| = y. 

Altitude of Maximum Concentration . 

The normalized altitude of maximum Ag is derived by 
maximizing Equation.. 5 with, respect to n . The resulting 
equation is transcendental, with no explicit expression for 
as a function of a and 8 . It can be written in the • 
following form: 

8 = (izHm) exp (2HEL) . 

1+n m a2 


( 6 ) 




Values of ri m that satisfy Equation 6 are plotted for various 
•values of 3 in Fig. 1. As expected, near the source (a -*• 0) 
surface bounce does not influence n m , which approaches one 
regardless of 6. This is to say Aq is a maximum at z = H. 

In the downwind direction 8 does strongly influence the 
altitude of maximum concentration. As ct -*■ n ra (1-8)/ (1+8) 
asymptotically. As a result, the. 8 *» 0 case, with no surface 
bounce, leaves = 1 for all a. For 6 = 1 or complete 

surface bounce., 0, which is to say the peak value descends 

to the earth's surface, a fact recently pointed out by 
Dumbauld (1976) .* Yet, n m only descends to the n = 0 plane 
if 8 = 1, and in that case n m = 0 for all a >_ 1. Finally, 
one notes n ra is a' strong function of 3 for a _> 1 , allowing 
it to serve as a primary detection characteristic of 6 according 
to Equation 6. 

‘ In application, H is estimated from standard plume rise 
formulas (Briggs, 1969), and the altitude, z, of maximum Aq 
is measured from the lidar vertical scans. Dividing z by H 
yields n ra . Next, cr z is estimated from the stability class 
and downwind distance x, using the graph on page 9 of Turner 
(1970) . Then a .is computed from Equation 3. Equation 6 is 
used to compute 8 from rim and a. A value of 8 is calculated 
for each a, and if the surface is simple, one expects them to 

t 

be approximately equal. A least squares method can be used 
to estimate 8 for the surface, provided the values are not 
excessively scattered. Measurements under various stability 
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conditions can examine the repeatability of- g for a given 
surface. 

Mean Altitude of Concentration . A second method of 
determining g is less sensitive and involves more 
computation than the above method, but it does offer an 
independent check on the model's internal consistency. 

The second method is to measure the. mean altitude (expec- 
tation value of the altitude) using the normalized 
concentrations, Aq, as a distribution function that weights 
n. The mean altitude, r^,- depends upon the sampling interval 
(from n ■ 0. to n * u) and is defined as: 


n„ = M i (ul • ' - (7) 

U M 0 (U) 


where ^(y) is-the n^ moment of n with respect to Aq, 
defined by 

%(y) -■= ^n^A d dn • (8) 


Performing the integration indicated by Equation 8 with * 
n - 0 and n - 1 gives the zeroeth and first order moments , 

M n (u> = /I [(1-8) erf (-i~) + erf(-HZi) + 6 erf (u+ljl , _ 

0 / 2 ■ SE SE /E 

and t • — 


(u) = fL 1 1+8) erf (-i-) + erft— ) - Serf 

/ * 1 /2a /2a J 


2 a 


2 

+a f(l+6)exp (-_!_) - exp(-iHZ^3 — ) - S exp )}, 

1 2a 2 2a 2 2a^ J 


( 10 ) 
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where erf(t) is the normal error function, defined 

erf(t) = / u exp (-u 2 ) du. (11) 

ST 0 

# 

Dividing Equation 10 by Equation 9 yields n^, as shown by 
Equation 7. The expression for n" u is rather formidable, 
but several limiting expressions derived from Equation 3 
and 10 make it more comprehensible. See Fig. 2. 

First, it is to be expected that near uhe source 
(a = 0) bounce plays an insignificant role, so that is 

independent of E, and furthermore, that ~ « 1 at a = 0. 

That is to say, the mean altitude of effluent at x = 0 is 
r) ss 1 or z = H. By noting that 

erf (t) * 1 (12) 

and taking the limit of Mq(u) and M^(y) as a 0 one 
discovers 

afj M 0 (y> = /2T , (13) 

Thus, as expected, their ratio, ^ = 1 at the source. In 
Fig. 2 all values of n u (and n m ) approach 1 as a approaches 0. 

Next, the behavior of the expression Mg and are 
simplified if u is allowed to approach infinity . Physically , 
this corresponds to the case of sampling the plume over an 
unlimited range of altitude from n = 0 to ri = ». In that 
case, the quantity n<» mfay be called the "true" mean altitude 
of the plume. The values of M Q and become 

JiS M 0 Cu) = M 0 (-> =yf [(l-BJerf^) + (1+6)] (14) 






and 


JiS M x (u) H [(l+6)erf (-ij + (1-sJ 

^ * 2 - /2ct 

+ a (1+6) exp (-_ ) . • (15) 

2a 2 

Their ratio, "n^, is plotted in Fig. 2 as the top pair of 
lines for the two cases 3-0 and B = 1. One sees that ri a 
expands without bound as a -*• 00 far downwind from the source. 
This occurs because Equation 14 , the denominator of Equation 7 , 
remains finite as a increases, but Equation 15, the numerator, 
has a final term that increases without limit as a does 
likewise. Thus, the true mean altitude of the normalized 
concentration rises above H the farther downwind the plume 
travels. 

Having examined the behavior of rf near the source (a =* 0) , 
and for the case of infinite vertical sampling, it remains . 
to consider the "measured" mean altitude of Ag , 7T . This 
measured value is achieved by allowing the measured Ag 
versus n profile to be truncated at an upper bound, n = y. 

True lidar systems will only scan to altitudes of n equal 
to a small number. Because all values of t\ larger than y 
are considered to contribute nothing to the computation of tT^, 
clearly 77^ < ri^ for all finite y. Use of Equations r 7 , - 9, and 
10 yields 77 . A specific typical case, y = 3, is plotted as 
the middle pair of lines in Fig. 2, illustrating the two 
cases B = 1 and 3 = 0. (The bottom pair of lines in Fig. 2 
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merely repeat for comparison the 6=0 and 8=1 cases of 
ri m previously plotted in Fig. 1) . 

Note that rT does not approach infinity as a increases, 
but instead remains finite throughout the entire range of a. 
Since only samples between n — Q and ri = ]i are considered., 
the value of must lie in the range 0 <. _<• u. Thus 

is always finite. 

Figure 2 reveals - an additional noteworthy fact. The 
measured mean altitude, is not nearly as strongly 

dependent upon 3 as n^, the altitude of maximum Aq / is. Thus 
as stated above; the measured mean value is not as sensitive 
a measure of 6 as n is. • 
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DOWNWIND FLUX 


The Pasquill-Gifford Model has no built-in pollutant 
removal mechanism. Thus, mass flux at the source Q, is 
equal to the mass flux through any vertical, semi-infinite 
(2 >0) plane normal to the wind at any distance downwind 
from the source. If 6 < 0, however, some material is 
removed from the plume and therefore downwind flux is 
affected. Define a normalized flux, 0(a), a function of 
normalized downwind distance, a, as:_ 

OQ OO 

0(a) = f f A dydn . - (16) 

y U — 00 

Substitution from Equation 3 shows 


0(a) = / / (£2)dydz 

VJ —OQ * 


(17) 


where C is a function of a implicitly. Returning to Equation 
16, and performing the y integration reveals 

/2tt o, r _ M 0 (~) 


0 (a) =• 


-Z S A dn = „ 

2ir a„ ^ 0 /2 tt 


(18) 


Since Mq («>) is /2n at a = 0, 0(a) = 1 at the source. But 
Mq (») is a function of 8 as well as a. Figure 3 reveals this 
dependence. 0(a) is plotted on the left ordinate, and it may 
be considered the fraction of material emitted at the source 
that is still airborne at a. The right ordinate plots 
Mq (°°) = 0(a). It is also useful in computing r^, for it 

is the normalizing factor. 



.0 


1.9 



Fig. 3. The normalized flux, 
effluent mass still airborne 


factor (M n (") ) . 



($(a)) or fraction of original 
at a, and the normalization 
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If $ = 0 the matter dispersed downward from H is 
.ultimately absorbed by the surface, but that borne upward 
remains aloft until other atmospheric processes invalidate 
the original assumptions of the model. Consequently, no 
more than half of the mass is ever removed from the plume 
by surface trapping. 


INFLUENCE Or 


ON SURFACE CONCENTRATION • 


Since 'humans live ' on the earth's surface, it is worth- 
while to investigate the significance of the - value being 
less than, one on surface concentrations. This is perhaps 
best seen in the .tore familiar basic model, rather than in 
the non-cimensional model. Thus going back to Equation 1, 
and setting z = C yields 

2 . 2 ' ' 


C = 


Q 


-«• * - — 


expf--^- --h— } r _+c) 
z ~v 


(19 ) 


'■Tz v - z 

For the Pascuiii-Gif ford Mode 1 5 = 1. Then the final (1-*-=) 
factor equals two. If the true value c f £ is less than cne , 
its influence on the surface concentration is clearly salutary 
but it only influences the concentration and not its geographic 
distribution. 



CONCLUSION 


In testing numerical models for plume dispersal with 
field experiments, it is common to reconcile discrepancies 
between predicted and measured values by adjusting source 
parameters. A simple adjustment in surface reflectance, 3, 
would often accomplish the same reconciliation without 
disturbing source variables one has no reason to change 
without knowledge of vertical wind profiles and surface 
roughness. Future lidar experiments may find it profitable 
to study sources that disperse over a variety of simple , yet 
diverse terrains, such as forest, grassland and water in 
order to investigate the usefulness of a parameter 8, and 
to see if it is itself independent of wind speed and stability 
class over a reasonable range of common conditions . If 
consistent values of 6 can be established, improvements 
in regional predictor models like the AQDM could be wrought. 
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Comparative Study of Various Text Editors and Formatting 
Systems 

ADA-029 050/ PAT 93 p. PC Not Avai!able/MF$3.00 

Explaining Energy: A Manual of Non-Style for the Energy 

Outsider Who Wants In 

LBL-4458/ PAT 78 p. PC$4.50/MF$3.00 

A Survey of Slate Legislation Relating to Solar Energy 
PB-258 235/ PAT 166 p. PC$6.75/MF$3.00 

Cost Estimating Handbook for Transfer, Shredding and 

Sanifary Landfilling of Solid Waste 

PB-256 444/ PAT 85 p. PC$5.00/MF$3.00 

Coal Liquefaction Design Practices Manual 
PB-257 541/ PAT 372 p. PC$1 0.50/M F$3.00 


HOW TO ORDER 


When you indicate the method of pay- 
ment, please note if a purchase order is not 
accompanied by payment, you will be billed 
an additional $5.00 ship and bill charge. And 
please include the card expiration date when 
using American Express. 

Normal delivery time takes three to five 
weeks. It is vital that you order by number 


or your order will be manually filled, insur- 
ing a delay. You can opt for airmail delivery 
for $2.00 North American continent; $3.00 
outside North American continent charge per 
item. Just check the Airmail Service box. If 
you’re really pressed for time, call the NTIS 
Rush Handling Service (703)557-4700. For a 
$10.00 charge per item, your order will be 
airmailed within 48 hours. Or, you can pick 
up your order in the Washington Informa- 
tion Center & Bookstore or at our Springfield 
Operations Center within 24 hours for a 
$6.00 per item charge. 


You may also place your order by tele- 
phone or if you have an NTIS Deposit Ac- 
count or an American Express card order 
through TELEX. The order desk number is 
(703) 557-4650 and the TELEX number is 
89-9405. 

Thank you for your interest in NTIS. Wc 
appreciate your order. 


METHOD OF PAYMENT 

n Charge my NTIS deposit account no. 

CH Purchase order no 

□ Check enclosed for $ 

0 Bill me. Add $5.00 per order and sign below. (Not avail- 
able outside North American continent.) 

□ Charge to my American Express Card account number 



NAME 

ADDRESS 

CITY. STATE. ZIP 


Card expiration date 

Signature 

Q Airmail Services requested 

Clip and mail to* 

NTS 

National Technical Information Service 
U.S. DEPARTMENT OF COMMERCE 
Springfield. Va. 2216! 

(703) 557-4650 TELEX 89-9405 




